Strong-coupling solution of the bosonic dynamical mean-field theory 
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We derive an approximate analytical solution of the self-consistency equations of the bosonic 
dynamical mean- field theory (B-DMFT) in the strong-coupling limit. The approach is based on a 
linked-cluster expansion in the hybridization function of normal bosons around the atomic limit. 
The solution is used to compute the phase diagram of the bosonic Hubbard model for different 
lattices. We compare our results with numerical solutions of the B-DMFT equations and numerically 
exact methods, respectively. The very good agreement with those numerical results demonstrates 
that our approach captures the essential physics of correlated bosons both in the Mott insulator 
and in the superfluid phase. Close to the transition into the superfluid phase the momentum 
distribution function at zero momentum is found to be strongly enhanced already in the normal 
phase. The linked-cluster expansion also allows us to compute dynamical properties such as the 
spectral function of bosons. The evolution of the spectral function across the transition from the 
normal to the superfluid phase is seen to be characteristically different for the interaction driven 
and density driven transition, respectively. 

PACS numbers: 71.10.Fd, 67.85.Hj 



I. INTRODUCTION 



Cold atoms in optical lattices provide a fascinating new 
class of interacting quantum many-particle systems J^ 
Due to the unprecedented precision of experimental tech- 
niques in this field it is now possible to simulate and 
experimentally test theoretical modelsi^"— In particular, 
experiments with bosonic atoms have revived the theo- 
retical interest in the properties of the bosonic Hubbard 
modeli^^— This model describes the quantum mechanical 
competition between the kinetic energy of lattice bosons, 
which is responsible for their Bose-Einstein condensation, 
and the repulsive interaction, which favors localization 
of the particles. The phase diagram of the bosonic Hub- 
bard model was first calculated by Fisher et al^ within 
a static mean-field theory derived from the atomic limit. 
With the formulation of the bosonic dynamical mean- 
field theory (B-DMFT ) ^^'^^ a comprehensive investiga- 
tion scheme for correlated lattice bosons in the thermo- 
dynamic limit has become available, which allows one to 
calculate also dynamical properties such as spectral func- 
tions of the interacting bosons. The B-DMFT is a ther- 
modynamically consistent, non-perturbative many-body 
approach which is applicable for all values of the input 
parameters, e.g., the interaction, density, and tempera- 
ture. It leads to a set of nonlinear equations which need 
to be solved self-consistcntly. An exact solution can be 
found only in special cases, e.g., for the Falicov-Kimball 
modeliiiii^ In general, the self-consistent equations have 
to be solved numerically or by employing approximate 
analytical methods. The experience with the fermionic 
DMFT— "— shows that both numerically exact (but com- 
putationally expensive) methods and approximate ana- 



lytical methods are important to gain insight into the so- 
lution of the complicated self-consistency equations. So 
far solutions of the B-DMFT equations had to be ob- 
tained fully numerically. Hu and Tongji^ and Hubener, 
Snoek, and Hofstetter— employed exact diagonalization 
(ED), and Anders et alr^^^ made use of continuous-time 
quantum Monte Carlo (CT-QMC) to solve the B-DMFT 
equations. Analytical or semi-analytical solutions of the 
B-DMFT equations did not exist up to now. 

In this paper we present an analytical strong-coupling 
solution of the B-DMFT derived by a linked-cluster ex- 
pansion (LCE)^i2i around the atomic limit. The method 
is analogous to the fermionic strong-coupling solver de- 
veloped by Dai, Haule, and Kotliar— for the fermionic 
DMFT. While in the fermionic case the strong-coupling 
expansion is unable to capture the low temperature Fermi 
liquid physics due to the existence of a characteristic low 
energy (Kondo) scale, there is no such limitation in the 
bosonic case. Our approach differs from previous strong- 
coupling expansions to the bosonic Hubbard model^^"— 
since they performed the expansion in the hopping am- 
plitude. 

The paper is organized as follows. We first introduce 
the B-DMFT and its self-consistency equations. Then 
we formulate the linked-cluster expansion and thereby 
derive a strong-coupling approximation to the B-DMFT 
equations. This is then applied to the Bethe lattice and 
the cubic lattice, both with coordination number z = 6, 
and to the Bethe lattice with z = oo. The phase dia- 
grams of the bosonic Hubbard, model calculated in this 
way are compared with those obtained from numerical 
solutions of the B-DMFT computed with EDii and CT- 
QMC fi^ii^ respectively, from numerically exact evalua- 



tions on a Bethc lattice?^ and from numerical results ob- 
tained by direct Monte Carlo simulations of the bosonic 
Hubbard model.— The momentum distribution functions 
and spectral functions of correlated lattice bosons in the 
normal and the Bose-Einstein condensed phase are also 
calculated. Finally we discuss possible extensions of the 
approach. 



II. CUMULANT EXPANSION IN THE 
BOSONIC DYNAMICAL MEAN-FIELD THEORY 

The B-DMFT is the bosonic counterpart to the well- 
established DMFT for lattice fermions described by the 
Hubbard model. Its derivation is described in detail in 
Ref. [ll|. Here we focus on a single species of bosons. 
The expansion presented below is easily generalized to 
the case of more than one type of boson. 

The bosonic Hubbard model is given by the Hamilto- 
nian 
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where h\ and hi are creation and annihilation operators, 
respectively, for a boson at a lattice site i, Uj is the hop- 
ping between lattice sites i and j, L'^ is the local interac- 
tion, and rii = b\hi is the number operator of the local 
occupation. In this paper we consider nearest-neighbor 
hopping, i.e., ty = —t for the nearest-neighbor sites i,j, 
and tij = otherwise. In the following we set the Boltz- 
mann constant fc^ and the lattice spacing a equal to 
unity. 



A. Local action of the B-DMFT 

In the B-DMFT the d-dimensional lattice problem ([1]) 
is replaced by an effective single-site ("impurity") prob- 
lem in which the local interaction U remains unchanged, 
but the rest of the lattice is replaced by two dynam- 
ical mean fields ("baths") corresponding to bosons in 
the normal state and in the Bose-Einstein condensate, 
respectively.— The time evolution of bosons on a partic- 
ular site z = is represented by the local Green function 



G(r) = -(T.b(r)bt(0)) 
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(2) 



where we used the imaginary time, finite temperature 
formalism and Nambu notation with 



(3) 



and the Bose-Einstein condensate (BEC) is described by 
the local order parameter 



= (Kt))s. 



(4) 



The impurity problem is defined by the local action 



Sloe = I dTb*{T){ — 
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+K / dT#^(T)b(T) + - dr dr' b^ (t)A(t - T')b(T'),(5) 

where fi is the chemical potential, k = X]i=^o ^«o is a lat- 
tice dependent parameter, and 
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is the condensate wave function, i.e., a dynamical mean 
field. The dynamical mean field corresponding to bosons 
in the normal state is represented by the hybridization 
function 



A(t-t') = 



An(r-r') Ai2(r - r') 
A2i(t-t') A22(r-r') 



(7) 



The dynamical mean fields $(t), Aii(t), and Ai2(t) are 
determined by the self-consistency equations 



^{t-t') = - ^ f,oio,(T.b,;(r)b](r'))(o) 

= J2t^otoM?ir-r') (8) 



and 



*=(b(r))(o). 
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Here the notation (■ • • )(o) indicates that the thermody- 
namic average is performed on a lattice with a cavity, i.e., 
with one site removed. We note that in equilibrium $(t) 
is constant. For finite dimensional lattices $ is related to 
the local BEC order parameter (|3]) by 
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The self-consistency loop is closed by introducing the 
self-energy in the Matsubara frequency representation 
through the k-integrated Dyson equation 



^(*^")^(, -^u;., + ^,)-^^'^^- 
and using the lattice Hilbert transform 
r^r \ [at I xTA^n + M-e 
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The latter equation links the local Green function to the 
self-energy for a specific lattice described by the non- 
interacting density of states Nq^s). The momentum de- 
pendent lattice Green function G(k, ia;„) is then given 
by 

-1 



G(k, iujn 



-iujn + n-ek 
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where ek is the dispersion relation of the non-interacting 
system and 5](iw„) is the self-consistent solution of equa- 
tions ([2l)-(fT2l). 

For a Bethc lattice with infinite connectivity (z = 
00)^1^ the self-consistency conditions reduce to the sim- 
ple expressions A(r — r') = t^G(r — r') and $ = {(/), (p*). 
In general, e.g., for a cubic lattice, the self-consistency 
equations (|T0]) - (fT2|) need to be solved numerically. 



B. Cumulant expansion 

In order to solve the impurity problem defined above 
we use the cumulant (linked-cluster) expansion in the 
dynamical mean fields An and A 12. The action ((S]) is 
further divided into two parts 



Here the superscript c indicates that only the connected 
part of the averages with respect to Sq is included. Now 
the partition function ([T7| can be calculated to the de- 
sired order in A. 

The above approximation is in the spirit of other 
strong-coupling expansions^^^ and becomes exact in the 
atomic limit (An — A12 = 0, (j> = 0). However, it should 
be stressed that it is not an expansion in the hopping 
amplitude but rather in the dynamical mean fields An 
and A 12. The fact that these fields are obtained self- 
consistently implies that all orders of the hopping ampli- 
tude contribute liii^S 

In the following we perform the cumulant expansion to 
second order in An and A12 in the partition function Z. 
Since the Green function is determined by the functional 
derivative 



"Jloc — Sq + S , 
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where 
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the diagonal element Gn (r — r') and off-diagonal element 
Gi2(t — t') are then of first order in An and A12: 



So = jM*{T)i— - f,)b{T) + -J^ dTUn{T){n{T) - 1) 

rfi 
K I dT$^(T)b(T) (15) 
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The partition function of the impurity problem Z is thus 
written as 



Furthermore, the local BEG order parameter is given by 



Z = Zo(e-^)o, 



(17) 
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where Zq is the partition function for the system de- 
scribed by So, and (...)o denotes the thermodynamic av- 
erage with respect to the action So- 

Now the exponential function appearing in the average 
is expanded, leading to an infinite series 

(e-^')o= 1-^ fdrfdr' (T,bt(r)A(r, r')b(r'))o 
^ JO Jo 

1 rP ffi ffS fP 

+ 4^/0'^^!/ "^^i/ 70"^^' (T.bt(ri)A(Ti,rO 



X b(T{)bt(r2)A(r2,r^)b(T^))o + 
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The series is then re-exponentiated with the help of cu- 
mulants (i.e., connected n-particle Green functions )^>^ 

(e-^')o= exp|-i^dr^dr' (T,bt(r)A(r, r')b(r'))S 

1 pP rP rP rP 
+ 4:2! y "^7 '^^V '^^^1 '^^'^ (r-bt(ri)A(ri, rO 
X b(r0bt(r2)A(r2,r^)b(r^))E; + ...}. (19) 
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(23) 
The thermodynamic averages are performed as 
(• • • )o ^ i^Tr{e-f^"° • • • ), with Ho = \Un{n-l)- ^in + 
K(<i>6^-f $*5). The trace is calculated over the eigenstates 
of Ho, which are obtained by an exact diagonalization of 
the Hamiltonian matrix which is represented in the oc- 
cupation number basis. Since the local Hilbcrt space of 
Ho for the bosonic impurity problem is infinite dimen- 
sional, the diagonalization has to be performed numeri- 
cally which, in principle, implies a further approximation. 
The Hilbert space has to be cut off in the occupation 
number of the impurity. The error introduced thereby 
can be controlled by performing calculations with differ- 
ent values of the cut-off and choosing the smallest cut-off 
value such that the results do not differ within the re- 
quired accuracy— 



III. APPLICATION OF THE LINKED-CLUSTER 
EXPANSION TO VARIOUS LATTICES 

In the following we apply the results of the LGE to 
the Bethe lattice and the cubic lattice, both with co- 
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. this work (B-DMFT with LCE) 
) numerically exact (cavity method), Ref. 28 
static mean-field theory, Ref. 8 




FIG. 1: Results for the Bethe lattice with coordination num- 
ber z = 6. Top: Dependence of the condensation tempera- 
ture Tbec on the interaction U at density (n) — 1. Bottom: 
Phase diagram n/U vs. t/U computed by different methods: 
B-DMFT solved with LCE (this work, T = O.lt), B-DMFT 
solved with ED (T = 0),— numerically exact evaluation (cav- 
ity method) (T < 0.25i),— and static mean-field solution 
(T = 0).— Inside the Mott lobes the system is Mott insu- 
lating and the occupation number is integer, while outside 
the system is superfluid. 



Snock, and Hofstetterjii and the static mean-field solu- 
tion of Fisher et al^ The static mean-field and the ED 
results were calculated at T = 0, whereas the results of 
the cavity method were obtained for T < 0.25t. The 
phase transition line ^/U vs. t/U only weakly depends 
on T at such low temperatures as can be seen in the up- 
per panel of Fig. [U where below TBEc{U)/t = 1 the curve 
is practically vertical. For this reason we conclude that 
the phase diagram presented in the lower panel of Fig. [1] 
is essentially the ground state phase diagram. 

The results shown in the lower panel of Fig. [T] demon- 
strate that the agreement between the two B-DMFT 
solutions is excellent. Namely, the blue circles (LCE, 
this work) are seen to lie practically on the red line 
(ED from Ref. [l7l ). Apparently the transition from the 
Mott-insulator to the superfluid is well described by the 
LCE approximation, which expands to first order in the 
dynamical mean field A(t). This is different from the 
case of the fermionic DMFT where the low temperature 
physics of the Hubbard model close to the metal- insulator 
transition can not be described by the strong-coupling 
approximation.— 

The value of the transition temperature Tbec(C^) ob- 
tained by the B-DMFT and the cavity method, respec- 
tively, is significantly lower than the results obtained by 
the static mean-field theory^ Since the B-DMFT cap- 
tures local fluctuations exactly we conclude that they are 
responsible for the lowering of Tbec{U) and the associ- 
ated increase of the size of the Mott lobes. 

For strong interactions the system is a Mott insula- 
tor for most values of the chemical potential /i. Upon 
lowering the interaction the system enters the superfluid 
phase with an order parameter cj) ^ 0. For the values of 
the chemical potential between the Mott lobes the super- 
fluid phase persists up to very large values of U. Since 
the LCE calculations were performed at a low but flnite 
temperature (T = O.lt), there is no superfluid phase at 
fi = U below t/U « 0.0001 (not discernible in the flgure). 



B. Cubic lattice 



ordination number z = 6, as well as to the Bethe lat- 
tice with infinite connectivity (z = (X)). Our results for 
the Bethe lattice with coordination number z < oo can 
benchmarked by the exact numerical solution based on 
the cavity method i^ 



A. Bethe lattice w^ith coordination number z = 6 

In Fig. [T] we show the results obtained with the LCE 
for the interaction dependence of the Bose-Einstein con- 
densation temperature Tbec{U), as well as for the phase 
diagram /z/C/ vs. t/U at T = O.lf. We also compare 
them with results from other methods: the exact numer- 
ical evaluation (cavity method) by Semerjian, Tarzia, and 
Zamponi^ the B-DMFT solution with ED by Hubener, 



1. Phase diagram 

The phase diagram of the Bose-Hubbard model for 
the cubic lattice obtained from the B-DMFT with the 
LCE and with CT-QMC, respectively, is presented in 
Fig. [21 These results are compared with the lattice quan- 
tum Monte Carlo (QMC) results,^ The LCE resuhs 
are shown for two different temperatures (T — 2t and 
T = 0.5t). It is evident that the size of the Mott lobes 
decreases with decreasing temperature. Upon lowering 
the temperature the computation of the phase boundary 
using the B-DMFT with the LCE was found to become 
more elaborate. As already noted in Ref. |l^ for the CT- 
QMC solver, the convergence of the DMFT cycle close 
to the phase transition is very slow and the initial guess 
of A(r) and $ has to be carefully chosen. 



this work (B-DMFT with LCE), T=2t 
this work (B-DMFT with LCE), T=0.5t 
B-DMFT with CT-QMC, T=0.5t, Ref. 18 
lattice QMC, T=0.025t, Ref. 29 
static mean-tield theory, Ref. 1 




FIG. 2: Phase diagram ^/U vs. t/U for the cubic lat- 
tice obtained from B-DMFT with LCE compared with the 
results obtained from B-DMFT with CT-QMC (data from 
Ref. \18^ . lattice QMC (data from Ref. l29l). and static mean- 
field theory.— 



Fig. [2] shows that there is a small quantitative dif- 
ference between the results obtained by different meth- 
ods. It is unlikely that these differences can be ex- 
plained by the different temperatures used in the compu- 
tations (the lattice QMC calculations^^ were performed 
at T = 0.025t, which is lower than the temperature used 
in the B-DMFT calculations). Indeed, at such low tem- 
peratures the temperature dependence of the phase di- 
agram is very weak, as discussed earlier for the z = 6 
Bethe lattice. Nevertheless, the overall agreement be- 
tween the results obtained from the three different meth- 
ods is clearly very good. As in the case of the z = 6 
Bethe lattice the local dynamical fluctuations described 
by the B-DMFT lead to an increase of the size of the 
Mott lobes compared to the static mean-field solution. 



2. Momentum distribution 

The momentum distribution function n(kx,ky,kz) of 
the normal phase, calculated at T = i, is found to have 
an interesting behavior close to the transition to the su- 
pcrfiuid phase. As shown in the upper panel and the 
inset of Fig. [3] the distribution n(fca.) = n{kx, 0,0) is 
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FIG. 3: Momentum distribution function for the cubic lattice 
in the normal phase calculated with the LCE at temperature 
T = t. Top: n(e) obtained for several values of the interaction 
close to the phase transition, which takes place at Uc = 28.5; 
inset: n{kx) for ky = kz = for the same parameters. Bot- 
tom: Comparison of the result for n(e) close to the transition 
{U = 1.13f/c) obtained by different methods. The results 
obtained with the 3'"'' order strong coupling expansion in the 
hop ping amplitude and the lattice QMC results are both from 
Ref. [23 



strongly peaked at fc^; = already in the normal phase. 
In the B-DMFT the momentum dependence of the mo- 
mentum distribution is expressed only through the non- 
interacting dispersion relation et (cf. Eq. ([T^l) ). There- 
fore, n(e) = n(ek) implicitly determines the momentum 
distribution. The plots in Fig.[3]show n(e) and n{kx) for 
different values of U upon approaching the phase transi- 
tion at constant density (n) = 1. 

The peak in the momentum distribution in the normal 
phase close to the phase transition was noted previously 
by Kato et al^ within QMC solution. The lower panel 
in Fig. [3] shows a comparison between n{e) obtained for 
the same parameters using different methods. As pointed 
out by Freericks et al^ the increase in the occupation 
at e = is an effect which is only partially described by 
a strong-coupling expansion in the hopping amplitude. 
The B-DMFT does capture this enhancement, and our 
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FIG. 4: Evolution of the spectral function (in arbitrary units) 
across the interaction driven phase transition at T = 2i; left 
column: phase transition at the tip of the Mott lobe, (n) = 1; 
right column: phase transition away from the tip, fi = 0.23C/. 
In both columns the bottom plot is for the normal phase, 
whereas the two upper plots are for the superfluid phase. The 
energy scale is plotted relative to the chemical potential which 
is at oj = 0. 



LCE results are in a very good agreement with the lattice 
QMC data of Rcf. H^. 



3. Spectral functions 

The B-DMFT approach with the LCE solver also al- 
lows one to investigate the behavior of the fc-integrated 
spectral function A{uj) = ——^^^ImGCkjUj) across the 
phase transition from the superfluid to the Mott phase 
(Figs. |4]and[5]). Since in our current implementation of 
the LCE the computations are performed on the imag- 
inary time or imaginary frequency axes, spectral func- 
tions at real frequencies have to be calculated by an- 
alytic continuation>2^ The spectral functions presented 
in Figs |4] and [5] were obtained by analytic continuation 
with Fade approximants. Calculations of bosonic spec- 
tral functions were also done with the functional renor- 
malization group^^i^l and in the variational cluster ap- 
proach (VCA)^i^ 

Here we focus on three distinct cases: (i) the interac- 
tion driven phase transition at the tip of the Mott lobe, 
keeping the ratio fi/U constant; (ii) the interaction driven 
transition at the bottom of the lobe, also with fi/U con- 
stant; and (iii) the density driven transition at constant 
interaction U. Due to the approximation introduced by 
the analytic continuation one can draw only qualitative 
conclusions about the spectral density in the region close 
to the chemical potential (e.g., one cannot reliably esti- 
mate the size of the gap). Nevertheless the qualitative 
behavior and the spectral weight transfer is well illus- 
trated and the difference between the three cases consid- 
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FIG. 5: Evolution of the spectral function (in arbitrary units) 
across the density driven phase transition at U = 3 If and 
T = 2i. The bottom plot is for the normal phase, whereas 
the two upper plots are for the superfluid phase. The chemical 
potential is at a; = 0. 



ered here is clearly visible. At the tip of the Mott lobe 
(case (i), left panel of Fig. U) an increase of the interac- 
tion leads to a symmetric shift of the spectral weight on 
both sides of the chemical potential. At the same time a 
Mott gap opens and two Hubbard bands are formed (see 
the bottom plot in the left panel of Fig. U) . The shape 
of the bands vaguely resembles the non-interacting den- 
sity of states iVo(e) for the cubic lattice. Away from the 
tip (case (ii), right panel of Fig. |4]) the shift of the spec- 
tral weight is not symmetric with respect to the chemical 
potential. The lower Hubbard band resides close to the 
chemical potential, whereas the upper Hubbard band is 
shifted to higher frequencies. A different behavior is ob- 
served in the density driven transition (case (iii). Fig. [5]). 
Upon increasing the chemical potential at constant in- 
teraction, the spectral function is shifted as a whole to 
lower frequencies, simultaneously forming a gap. 



C. The Bethe lattice with z = oo 

The phase diagram for the z = oo Bethe lattice is pre- 
sented in Fig. [6l At sufficiently high temperatures (e.g., 
T = 0.6t as in Fig. ^ the LCE gives convergent results 
both for the superfluid and normal phases near the phase 
transition. However, at temperatures below 0.4t we have 
not been able to find a convergent solution in the super- 
fluid phase around the tip of the second Mott lobe. The 
iterations converge either to (p — (normal phase), or 
to a solution with (f) ^ but with a non-concave, and 
hence unphysical,— G'ii(r). The results for T = OM are 
shown in Fig. [SI where the solution at fi/U around the 
first Mott lobe converges both in the normal and the su- 
perfluid phase, thus making it possible to calculate the 




FIG. 6: Phase diagram ^/U vs. t/U for the z = oo Bethe 
lattice at two different temperatures. At T — 0.3i it was not 
possible to determine the phase boundary of the second lobe 
in the region 1.26C/ < fi < 1.8U (see the discussion in the 
text). 



phase boundary. In the range 1.26U < fi < 1.8C/ a con- 
vergent solution was only obtained in the normal phase 
(0 = 0), i.e., it was not possible to determine the phase 
boundary of the second lobe completely. As the temper- 
ature is lowered, the range of the chemical potentials for 
which we did not obtain a supcrfluid solution increases. 
For temperatures below 0.3t we did not even obtain so- 
lutions with non-zero superfluid order parameter around 
the tip of the first Mott lobe. Upon further lowering the 
temperature, the region of convergence of the method in 
the superfluid phase is reduced to the values of ^ near the 
edges of the lobes. At this moment it is not clear whether 
the absence of a solution in the superfluid phase in the 
z = oo Bcthe lattice for some chemical potentials at low 
temperatures is a consequence of the strong-coupling ap- 
proximation to the B-DMFT, or the B-DMFT itself. This 
is an open question which needs to be answered in the 
future. Such problems did not occur for the other lattices 
investigated here. 



IV. SUMMARY 

We developed an analytical approximation scheme to 
solve the B-DMFT equations for correlated lattice bosons 



in the strong-coupling limit. The solution makes use of a 
linked-cluster expansion in the hybridization function of 
normal bosons around the atomic limit. Explicit results 
were obtained for the Bose-Hubbard model on the cubic 
lattice and the Bethe lattice with connectivity z = 6 and 
z = oo, respectively. Remarkably good agreement with 
numerical solutions of the B-DMFT equations obtained 
with exact diagonalizationii, continuous-time quantum 
Monte Carle — , and direct lattice QMC calculations^^ 
was found. This agreement demonstrates that the strong- 
coupling solution derived here provides a correct descrip- 
tion of the physics of correlated bosons. The method is 
computationally inexpensive and, with a good choice of 
the initial guess of the parameters, usually leads to a fast 
convergence of the iteration of the self-consistency equa- 
tions. The Bethe lattice with infinite connectivity is an 
exception which still requires further investigation. 

We also employed the linked-cluster expansion to cal- 
culate the momentum distribution function of normal 
bosons close to the phase transition as well as the bosonic 
spectral function in the normal and superfluid phase. 

The approximation scheme presented in this paper can, 
in principle, be systematically improved by the inclusion 
of higher order terms. However, the non-interacting limit 
can only be reached if terms up to infinite order are in- 
cluded, e.g., by an appropriate resummation. This has 
been achieved for fermions by the non-crossing approxi- 
mation (NCA)4i The fundamental problem of the NCA, 
namely its failure to describe the low temperature Fermi 
liquid regime adequately owing to the existence of a char- 
acteristic coherence scale (the Kondo temperature) , may 
be absent in the case of bosons where such a coherence 
scale does not exist. For that reason it should be clar- 
ified whether it is possible to construct a renormalized 
expansion for correlated bosons which is applicable for 
all temperatures and interaction strengths. 
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